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Abstract: 

We have constructed a Monte Carlo generator^ for lowest-order predictions for the 
processes 77 — > 4/ and 77 — > 4/7 in the Standard Model and extensions thereof by an 
effective 77H coupling as well as anomalous triple and quartic gauge-boson couplings. 
Polarization is fully supported, and a realistic photon beam spectrum can be taken into 
account. For the processes 77 — > 4f all helicity amplitudes are explicitly given in a 
compact form. The presented numerical results contain, in particular, a survey of cross 
sections for representative final states and their comparison to results obtained with the 
program package Whizard/Madgraph. The impact of a realistic beam spectrum on 
cross sections and distributions is illustrated. Moreover, the size of various contributions to 
cross sections, such as from weak charged- or neutral-current, or from strong interactions, 
is analyzed. Particular attention is paid to W-pair production channels 77 — > WW — > 
4/(7) where we investigate the impact of background diagrams, possible definitions of the 
W-pair signal, and the issue of gauge-invariance violation caused by finite gauge-boson 
widths. Finally, the effects of triple and quartic anomalous gauge-boson couplings on 
cross sections as well as the possibility to constrain these anomalous couplings at future 
77 colliders are discussed. 
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^The corresponding Fortran code can be obtained from the authors upon request. 



1 Introduction 

A photon (or 77) collider [1] as an option at a future e + e~ linear collider extends the 
physics potential of such a machine substantially. It provides us with information about 
new physics phenomena, such as properties of the Higgs boson or of new particles, which 
is in many respects complementary in the e + e~ and 77 modes (see, e.g., Refs. [1,2] and 
references therein). For instance, Higgs-boson production in the s-channel [1-4] becomes 
accessible which is highly suppressed in e + e~ annihilation, but very interesting for Higgs 
bosons of relatively large mass. Owing to its large cross section, W-pair production is of 
particular interest at a 77 collider. For instance, it can be used for precision tests of the 
gauge sector of the Electroweak Standard Model (SM). While the reaction e + e~ — > WW 
depends on the gauge-boson couplings ZWW and 7WW, the corresponding reaction at a 
77 collider, 77 — > WW, is sensitive to the gauge-boson couplings 7WW and 77WW. At 
an e + e~ collider the coupling 77WW is only directly accessible through the bremsstrah- 
lung process e + e~ — > WW7 which is suppressed by a factor a(0)/7r w.r.t. the non-radiative 
process e + e~ — > WW. Therefore, the sensitivity to the anomalous 77WW coupling in 
the 77 mode is expected to be an order of magnitude better than in the e + e~ mode. The 
precision for the measurement of the 7WW coupling is comparable in both modes [5]. 

Since W bosons decay into fermion-anti-fermion pairs, the actually observed final 
states of 77 — > WW involve four fermions. Lowest-order predictions for 77 4/ pro- 
cesses (with monochromatic photon beams and leptonic or semileptonic final states) were 
discussed in Refs. [6,7]. With the exception of Ref. [8], the existing analyses on non- 
standard couplings at a 77 collider, which focus on anomalous triple gauge-boson cou- 
plings (ATGC) [5,9], on anomalous quartic gauge-boson couplings (AQGC) [3,10], on 
CP- violating gauge-boson couplings [11], and on effects of strongly interacting longitudi- 
nal W bosons [12], treat W bosons as stable. In the above studies radiative corrections 
were not properly taken into account either. 

In view of the achievable experimental precision, the theoretical predictions should 
include, however, both radiative corrections and the decays of the W bosons, i.e., the full 
process 77 — > 4/ should be considered. The lowest-order amplitudes to the processes 
77 — > 4/ include diagrams with two resonant W bosons ("signal diagrams") as well as 
diagrams with only one or no W resonance ("background diagrams"). Compared to the 
doubly-resonant diagrams, such singly-resonant and non-resonant diagrams are suppressed 
by roughly a factor Tw/Mw and (Ty? / Mw) 2 , respectively. On the other hand, radiative 
corrections modify the theoretical predictions by contributions of 0(a). Since r w /M w 
is also of 0(a), both the off-shell effects and the radiative corrections contribute to the 
77 — > WW — > 4/ cross section in the same perturbative order. Hence, it is a promising 
approach to calculate the lowest-order matrix elements to the full 77 — > 4/ process and 
take into account only the doubly-resonant part of the 0(a) corrections to the reaction 
77 — > WW — > 4/. This approach is known as double-pole approximation (DPA) where 
only the leading terms in the expansion of the amplitude around the poles of the resonant 
W propagators are included in the calculation. The overall accuracy is then 0(a) relative 
to the lowest-order prediction for the doubly-resonant W-pair signal. This approach has 
been described in Ref. [13] in some detail and successfully applied to W-pair production 
in e + e~ annihilation in different versions [14-18]. 
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The doubly-resonant virtual corrections consist of factorizable and non-factorizable 
contributions. The former comprise the corrections to on-shell W-pair production [19- 
21] 1 and the decay of on-shell W bosons [23], the latter account for soft-photon exchange 
between the production and decay subprocesses [24]. Although the basic building blocks 
for the virtual corrections exist in the literature, the combination into a complete set 
of 0(a) corrections in DPA has not been done yet. The evaluation of bremsstrahlung 
processes 77 — * 4/ + 7 considered in this work serves as a further building block for the 
0(a) corrections in DPA. 

The calculation of lowest-order predictions for the processes 77 — > 4/ and 77 — > 4/7 
described in this paper closely follows the approach of Ref. [16] for the e + e~ case. Fermion 
masses are neglected everywhere assuming that all mass singularities are avoided by phase- 
space cuts. Following Ref. [16] we construct an event generator based on multi-channel 
Monte Carlo integration [25] using adaptive self-optimization [26]. A realistic photon 
spectrum can optionally be included in the parametrization of CompAZ [27] which is 
based on the results of Ref. [28]. In addition to the pure SM case, we calculate also 
the lowest-order matrix elements including anomalous couplings. Specifically, we con- 
sider ATGC [29,30] corresponding to SU(2) xU(l)-gauge-invariant dimension-6 operators 
and genuine AQGC [31,32] which respect electromagnetic gauge invariance and custodial 
SU(2) C symmetry. Finally, we include an effective 77H coupling in our analysis in order 
to study the Higgs resonance in the reaction 77 — > H — > WW/ZZ — * Af. 

The performance of the Monte Carlo generator is illustrated in a detailed numerical 
discussion. We present a survey of cross sections for a set of representative 4/ and 4/7 
final states and compare them with results obtained by the combination of the programs 
Whizard [33] version 1.28 and Madgraph [34]. We find agreement within statistical un- 
certainties. Moreover, we discuss some invariant-mass, energy, and angular distributions. 
We illustrate the impact of the realistic beam spectrum on cross sections and distribu- 
tions, and study the size of various contributions to cross sections, such as from weak 
charged- or neutral-current, or from strong interactions. In view of the aimed precision 
calculation for 77 — > WW — > 4/(7), which will include radiative corrections in DPA, we 
elaborate on the possibility to define a proper W-pair production signal. In this context, 
we investigate a naive signal definition which is based on doubly-resonant diagrams only 
and thus violates gauge invariance. We find that this naive definition, which works well 
for the e + e~ case (cf. discussion of the so-called "CC03" cross section in Ref. [35]), is 
not satisfactory for 77 — * WW — > 4/, but requires a proper gauge-invariant definition 
of the residue on the double resonance. For some 77 — > 4/ and 77 — > 4/7 processes 
we investigate the effects of gauge-invariance violation by introducing gauge-boson decay 
widths. To this end, we compare results obtained with the naive introduction of constant 
or running widths with the gauge-invariant result obtained with the so-called "complex- 
mass scheme" [16]. Finally, we discuss the effects of ATGC and AQGC on cross sections. 
In particular, we estimate the bounds on anomalous couplings that could be set by a 77 
collider, by assuming a photon spectrum and luminosities expected for the 77 option at 
TESLA. 

The paper is organized as follows. In Section 2 we describe the calculation of the 
helicity amplitudes in the SM and give explicit results for 77 — > 4f. In Sections 3 and 4 

1 Radiative corrections to on-shell W-pair production, 77 — ► WW, were also considered in Ref. [22]. 
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we describe the calculation of the amplitudes with anomalous gauge-boson couplings and 
an effective 77H interaction. Section 5 provides some details about the numerical Monte 
Carlo integration. In Section 6 we present the numerical results as outlined above. Our 
summary is given in Section 7. 

2 Analytical results for amplitudes in the Standard Model 

2.1 Notation and conventions 

We consider reactions of the types 

7(fci, Ai) + 7O2, A 2 ) -> /i(pi, 0-1) + / 2 (p2, cr 2 ) + / 3 (ps, 0-3) + /4(P4, 0-4), (2.1) 

7(fci, Ai) + 7O2, A 2 ) -> /i(pi, (Tl) + / 2 (P2, ^2) + /3(P3, ^3) + /4(P4, 0-4) + t(P5, 0-5) ■ (2.2) 

The arguments label the momenta fcj, and helicities A&, 0"^ (which take the values ±1/2 
in the case of fermions and ±1 in the case of photons) of the corresponding particles. 
We often use only the signs to denote the helicities. The fermion masses are neglected 
everywhere. 

For the Feynman rules we follow the conventions of Ref. [36]. We extend the usual 
linear gauge-fixing term in the 't Hooft-Feynman gauge by a non-linear part according 
to Refs. [16,20,37] such that the vertex jWcf) vanishes, where are the would-be Gold- 
stone bosons corresponding to the W bosons. Note that this also affects the gauge-boson 
couplings 77WW and 7WW. The corresponding Feynman rules relevant for 77 — > 4/(7) 
in lowest order can be found in Ref. [16]. Since we neglect fermion masses, the would-be 
Goldstone bosons do not couple to fermions and do not occur in the Feynman graphs of 
the SM amplitudes to 77 — > 4/(7), which leads to a considerable reduction of the number 
of Feynman diagrams. 

2.2 Classification of final states for 77 — > 4/(7) 

The final states for 77 — > 4/ and 77 — * 4/7 can be classified similarly to the processes 
e + e~ — > 4/ and e + e~ — > 4/7 [16]. In the following, / and F are different fermions (/ 7^ F), 
and /' and F' denote their weak-isospin partners, respectively. We distinguish between 
states that are produced via charged-current (CC), via neutral-current (NC) interactions, 
or via both interaction types: 

(i) CC reactions: 

77 - ff'FF' 

(ii) NC reactions: 

(a) n^ffFF (NC40 family), 

(b) 77 - //// (NC2-40 family), 

(iii) Mixed CC/NC reactions: 

77 -»• ///'/' (mix71 family). 



(CC31 family), 
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The radiation of an additional photon does not change this classification. Following 
Ref. [38] we give the names of the process families in parentheses where the numbers 
correspond to the number of Feynman diagrams involved in unitary or non-linear gauge 
(for processes without neutrinos in the final state, not counting gluon-exchange diagrams). 

Since the matrix elements depend on the colour structure of the final state we further 
distinguish between leptonic, semi-leptonic, and hadronic final states. Keeping in mind 
that we neglect fermion masses, omitting four-neutrino final states, and suppressing re- 
actions that are equivalent by CP symmetry we end up with 17 different representative 
processes which we have listed in Table 1. 



nnai state 


leaction type 


77 — >■ 


leptonic 








iNL^aj 








e e~*~/i fj,^ 




IN O ^ D ) 


e e e e 




CC/NC 


e~e + z/ e ^e 


semi-leptonic 


CC(c) 


e~z/ e ud 




NC(a) 


u c u e \x\x 






v c v e dd 






e~e + uu 






e~e + dd 


hadronic 


cc 


udsc 




NC(a) 


uucc 




NC(a) 


uuss 




NC(a) 


ddss 




NC(b) 


uuuu 




NC(b) 


dddd 




CC/NC 


uudd 



Table 1: Set of representative processes for 77 — * 4/(7). 

Since the photons are polarized after Compton backscattering, final states that are 
flavour equivalent up to a CP transformation need not necessarily yield the same cross 
section if the convolution over a realistic photon beam spectrum is included. However, 
as we neglect fermion masses, this is only relevant for the semileptonic CC processes 
77 — > e~z/ e ud(7) and 77 — > z/ e e + du(7). 

2.3 Lowest-order amplitudes for 77 — > 4f 

2. 3. 1 Construction of matrix elements 

The amplitudes for the processes 77 — > 4/ are constructed by attaching the two 
incoming photons in all possible ways to the corresponding diagrams with four external 
fermions as shown in Figure 1. The matrix element of the generic diagram in Figure 1, 
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Figure 1: Generic diagram for the process 77 — > 4/ where the photons 71,72 couple to 
the fermions fi, . . . , / 4 and the gauge boson in all possible ways. 

where two fermion lines are linked by a gauge boson V, can be written as 

■M-xlxly 4 (^i,Pj,Qj) = 4e 6 ai - a2 5 a3i ^ a4 g v ^^gy^^A x 1 iX i 2 V (ki,pj,Qj), (2.3) 

where hi, Pj and Qj (i = 1,2; j = 1,..,4) stand for the momenta and relative electric 
charges of the particles, respectively. The coupling factors (relative to the electric unit 
charge e) read 

s I 3 ■ 

a a - = -O- a a - =-—Q + w>% 5 

y-yfifi y Zfifi r r q 



Sw/^ - ^ ^ ~ e' (2 - 4) 



where 1^ = ±1/2 denotes the weak isospin of the fermion fi and g s the strong coupling 
constant. The weak mixing angle is defined by 

cos6> w = c w = — ^, s w = t/l - cl. (2.5) 



Quark mixing is neglected everywhere, i.e. we set the CKM matrix equal to the unit 
matrix. The auxiliary functions v are calculated within the Weyl-van-der-Waerden 
(WvdW) formalism following the conventions of Ref. [39]. The WvdW spinor products 
are defined by 

(pq) = e AB PA q B = 2^/p^ (e'^ cos ^ sin ^ - e"^ cos ^ sin ^ , (2.6) 

where pa, Ha are the associated momentum spinors for the momenta 

P^ = Po (1) srn @p cos 4>pi srn s i n cos ^p) 5 

<? M = 9o(l, sin 9 cos0 g , sin 9 sin g , cos# g ). (2.7) 
Moreover, we define the shorthands 

(piPkPj) = P it A P k B Pj,B = Pi,APkPkPj,B = (PiPk)*(PjPk), 
(pilPt + P m ] Pj ) = (piPiPj) + (piPmP,), (2.8) 
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where Pk,i, m are light-like momenta, i.e., p\ = pf = = 0. In the following, the 
denominators of the gauge-boson propagators are abbreviated by 



Pvip) 



y = 7 ,Z,W,g, M- M g 0. 



p 2 -Mf 

The introduction of the finite width is described in Section 2.4 below. 
The auxiliary functions A^^ y explicitly read 

A-- y (k i ,p j ,Q j ) = ((pm)*f 



(2.9) 



x \-Q\ 



(piP 2 )*(p3P4)Pv(P3+Pi) 



-Q1Q3 



(Pi +P2- k 1 ) 2 Py(p 1 +p 2 ~ h) 

(km)*(k\P 2 )*( k 2 p 3 ) * ( k 2 p 4 ) * 



(km)* (hp 2 )* (km)* (km)' 
+ Qz{Qi-Q 2 )Pv{pi+P2) 

-(P2P4)*(PiP 2 ) + (km) * (km) M V P V ( Pl +p 2 - ki) 



X 



( km ) * ( km ) * ( km ) * ( km ) * 



+ (ki <-> k 2 



+ (Qi - Q2?Pv{pi + p 2 )Pvipz + Pi) 



■(P2Pi) 



(pm)* (pm) (p3Pa) + M v (pm) 

2 (km)* (km)* ( hp 2 ) * ( k 2 p 4 ) * 



+ M v P v (p 1 +P2-k 1 ) 



(km) (km) 



(km)* (km)* 

+ ({Pl, <5i;P2, Q2} ^{P3, Qs; P4, Q4}) |, 



+ QlPviPs + Pa) 



(pm)*(km)(k 2 [P 2 + P4]p 3 ) 

(km)*(k 2 pi)(p 2 +P3 +P4) 2 

(k 2P2 )*(pm)m[Pi + Ps\ki) 



+ QlQ 2 P V (P3+P4) 



(km)*(k 2 p 2 )(pi +P3+P4) 2 
(p 2 [K 1 -P 1 }p 3 )(p 4 [K 1 -P 3 ]p 1 ) 



(km)* (km)* (km) ( km ) 



+ (Q 2 -Qi)Pvm + pd 



(pm)*(pm) 

(km)*(hpi) 



2 (km)* (km) 1 (km)* (km) 



1 



f |(<?2 - Qi) 2 ^y(pi + P2)iV(p 3 + p 4 ) 



x 



(P2P4)*(pm)((p2[K 1 - Pi]p3)(p±[Ki - p 3 ]pi) - M v (pmY(pm)) 



( km )*(km)*(km)( km ) 
+ [-Qi + (Qi - Q 2 )2(k lPl )Pv(Pi+P2)} [Qa + (Q3 - Q4)2(km)P v ( P3 +p 4 )] 

((p 2 [K l -P l }p 3 )) 2 Py( Pl +p 2 -k 1 ) 



X 



(km)*(km)*( km ) ( km) 

+ ({Pl,Ql,P2,Q 2 }^{P3,Q3;P4,Q4})- 

The other auxiliary functions A^J^y follow from the relations 

A \^\2,v(ki-,Pj-,Qj) — A x 1 iX i 2 y(ki,pj,Qj 
A\\\ 2 y(ki,Pj,Qj) — A x ri i a x i 2 y(k i ,pj,Qj) 



(2.10) 



{Pl,Ql}^{P2,-Q'2} ' 

{pz,Qz}^{pi-Qi\ 



(2.11) 



and 

A-xl'-xlyikiiPj, Qj) = ^-xlxl,v(h,Pj,Qj) ■ (2-12) 

The last relation expresses a parity transformation. Note that the operation of com- 
plex conjugation in Eq. (2.12) must not affect the gauge-boson widths in the propagator 
functions Py which will be introduced in Section 2.4. 

The calculation of the helicity amplitudes for 77 — > 4/7 proceeds along the same lines. 
The result, however, is quite lengthy so that we do not write it down explicitly. 



2.3.2 Squared amplitudes for leptonic and semi-leptonic final states 



The result for leptonic and semi-leptonic final states follows immediately from the 
generic amplitude (2.3). The gauge boson cannot be a gluon in this case, and the sum 
over the colour degrees of freedom in the squared matrix elements trivially leads to the 
global factors N£ t = 1 and iV s c emilept = 3. Note that for NC diagrams the result for the 
amplitude is much simpler than for CC diagrams, since all terms in Eq. (2.10) involving 
a factor (Q 1 — Q 2 ) or (Q 3 — Q 4 ) drop out. Most of these terms originate from diagrams 
where a photon couples to a virtual W boson. 

The explicit results for the colour-summed squared matrix elements read 



E \Mcc\ 

colour 

E l-^NC(a)| 
colour 

E \M NC (b)\ 

colour 

E 1-McC/Ncl 
colour 

where we use the shorthand 



= N' 



A^IMwl 2 , 
N C \M NC \' 
M NC - 
M NC - 



[■^ NC ]{pi,Ql,<Tl}~{ps,Q3,»S} 



(2.13) 
(2.14) 

(2.15) 

(2.16) 

(2.17) 



M^c = E M v 

y=7,z 

and suppress the helicity indices and the dependence on momenta and relative charges. 
The relative signs account for interchanging external fermion lines. 



2. 3. 3 Squared amplitudes for hadronic final states 

Next we consider purely hadronic final states, i.e., the cases where all final-state 
fermions are quarks. This renders the summation of the squared matrix elements over the 
colour degrees of freedom non-trivial, and in addition gluon-exchange diagrams appear. 
Since gluon-exchange diagrams require two quark-anti-quark pairs in the final state they 
do not appear in CC processes. For CC processes there is only one possibility for the 
colour flow, and the summation over the colour degrees of freedom leads to an overall 
factor iV£ adCC = 3 2 = 9 to the squared matrix elements as given in Eq. (2.13). 

For NC reactions we have to compute the sum of pure electroweak (ew) and gluon- 
exchange (QCD) matrix elements, 

\AC\C2CZCA \ AClC2C3C4 I A / fClC2C3C4 / r> -I q\ 
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where q denotes the colour indices of the quarks. The electroweak diagrams are diagonal 
in colour space and therefore read 



KA C1C2C3C4 
• /vl NC(a),had,cw 



-^NC^ClC2^C 3 C4! 



^NC(b)!had,ew - MncS Ci c 2 S C3C4 - [Mnc] {puQi 



(2.19) 



The gluon-exchange diagrams are obtained from the generic formula (2.3) by inserting 
the corresponding generators, A a /2, of the gauge group SU(3), 



\A C1C2C3C4 

- /vl NC(a),had,QCD 
• /vl NC(b),had,QCD 



= M g \x a clC2 \ a C3C4 - [M 



s J {pi,Qi^i}-{P3:Q3,<T3} ^c 3 c 2 K lC4 - (2.20) 



The matrix element M. g is defined by Eq. (2.3) with V — g. 

Carrying out the colour sum using the completeness relation for the Gell-Mann ma- 
trices, 



yields 

,had I 

colour 



9\M 



NC 



\ a \a 

A ij A kl 



+ 2\M 



(2.21) 



gi ' 



E l-^NC(b),had| 2 = 9|^ NC | 2 + 9|[^ N c] {piiQli(Tl}M{p3iQ3i(T3} | +2\M, 



colour 



+ 2 
4 
3 



[M 



S-l {Pl,Ql,CTl}<->{p3,Q3,0-3} 



+ - Re { M g M* g 



8Re{M g [M* NC ] {puQl 



6Re{M NC [M*v C ] {puQuai} „ M3ia3} } 
^ -8Re ^M nc [m; 

CTl } +^{P3,Q3,o-3}} ' 



{pi,Ql,0l} ^{p 3 ,Q3 



{Pl,Ql,o-l}^{p 3 ,Q3 



^3}} 

(2.22) 



All squared matrix elements of this section have been numerically compared with 
results obtained with the program Madgraph [34] at several phase-space points, and 
perfect agreement has been found. 

2.4 Implementation of finite gauge-boson widths 

We have implemented the finite widths of the W- and Z-boson propagators 2 in four 
different ways: 



• fixed width in all propagators: 

Pv(p) = 



/r' - My ■ \My\y' 



(2.23) 



2 We have also supplemented the explicit gauge-boson masses appearing in the numerators of Eq. (2.10) 
by the corresponding widths, because these mass terms originate from denominators upon combining 
different diagrams. 
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step width (fixed width in time-like propagators): 



= ^-j^+Uw) - (2 ' 24) 

• running width in time-like propagators: 

PyW = v 2 -Ml + if(T v /M v )6^y (2 ' 25) 

• complex-mass scheme [16]: complex gauge-boson masses are used everywhere, i.e. 



My — iMyTy instead of My in all propagators and couplings. This results in a 
constant width in all propagators, 

and in a complex weak mixing angle 



Cw_i Sw " M|-iM z r z • 10 



The virtues and drawbacks of the first three schemes are discussed in Ref. [40]. All 
but the complex-mass scheme, in general, violate SU(2) gauge invariance, the step- and 
the running-width schemes also violate electromagnetic U(l) em gauge invariance, which is 
preserved by using a fixed width. As known from many examples in e + e~ physics [16,40, 
41], gauge-invariance-violating effects, in particular when enhanced by factors p 2 /M v as 
in the running-width scheme, can lead to totally wrong results. Furthermore, the violation 
of U(l) em gauge invariance also causes a dependence of matrix elements and cross sections 
on the gauge chosen for external photons. In e + e~ — > 4/ and e + e~ — > 6f this problem 
does not occur since no external photons are involved. 

The complex-mass scheme, which was introduced in Ref. [16] for tree-level calculations, 
preserves gauge invariance and thus all Ward identities which rule gauge cancellations. Its 
application is particularly simple for 77 — > 4/(7) in the non-linear gauge. In this case, no 
couplings involving explicit gauge-boson masses appear, and it is sufficient to introduce 
the finite gauge-boson widths in the propagators [cf. Eq. (2.26)] and to introduce the 
complex weak mixing angle (2.27) in the couplings. 

For CC processes 77 — > 4/(7) with massless fermions, the fixed- width (FW) approach 
in the non-linear gauge and the complex-mass scheme (CMS) are practically equivalent, 
because all Feynman graphs are proportional to e 4 /s^, (e 5 /s^,) and gauge-boson masses 
appear only in propagator denominators. In this case the corresponding amplitudes in 
the two schemes differ only in the global factor 4,fw/4,cmsi where s Wj fw an d s w ,cms 
are the values of s w in the different schemes, i.e., s Wj fw is derived from the ratio of real 
gauge-boson masses and s Wj cms from complex masses. Thus, both squared amplitudes 
are gauge invariant and are equal up to the factor |s w ,fw/ s w,cms| 4 which is equal to 1 up 
to terms of 0{T^/M^). 
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For NC and CC/NC processes a similar reasoning can be used to show that the fixed- 
width approach does not violate gauge invariance in 77 — > 4/(7) for massless fermions. 
The trick is to apply the above argument to gauge-invariant subsets of diagrams. For NC 
diagrams with photon exchange, which is the (gauge-invariant) QED subset of diagrams 
(Figure 1 with V — 7), there is nothing to show. The sum of NC diagrams of type NC(a) 
with Z-boson exchange (Figure 1 with V = Z) again involves c w and s w only in a global 
coupling factor (per helicity channel); the remaining dependence on the gauge-boson 
masses is located in the propagator denominators. Thus, the subamplitudes of the fixed- 
width and the complex-mass scheme are again identical up to a global factor and both 
preserve gauge invariance and Ward identities. For NC processes of type NC(b) a second 
class of diagrams exists (Figure 1 with V — 7, Z and external fermions interchanged). 
This new class of diagrams forms a gauge-invariant subset because of the different flow of 
fermion numbers. Thus, the reasoning for type NC(a) applies to both classes of diagrams 
of NC(b) reactions. The same argument is also valid for the subset of CC diagrams in 
mixed CC/NC reactions. 

In summary, we have argued that the use of naive fixed gauge-boson widths does not 
lead to gauge- invariance violations in amplitudes for 77 — > 4/(7) as long as fermions are 
massless and the non-linear gauge with vanishing 7W0 coupling (or the complex W-boson 
mass in this coupling if the 't Hooft-Feynman gauge is chosen) is used. The corresponding 
squared amplitudes agree with the ones of the (gauge-invariant) complex-mass scheme up 
to terms of 0(r w /M w ), for CC processes even up to terms of 0(Y'^ N /M^ r ). 

2.5 W-pair signal diagrams and double-pole approximation 

The diagrams to CC and CC/NC processes comprise graphs with two, one, or no 
internal W-boson lines that can become resonant, similar to the situation for e + e~ — > 
WW — > 4/ (see Refs. [35,42] and references therein). It is interesting to investigate the 
possibility to define an amplitude for the W-pair signal based on doubly-resonant con- 
tributions only, because such an amplitude is much simpler than the full amplitudes for 
four-fermion production and is universal (up to colour factors) for all relevant 4/ final 
states. Moreover, this study is an important exercise for the calculation of radiative correc- 
tions to 77 — > WW — > 4/ in the so-called double-pole approximation (DPA), where only 
doubly- resonant contributions are taken into account. Taking simply all doubly- resonant 
diagrams, of course, yields a result that is not gauge invariant. Nevertheless in the e + e~ 
case the lowest-order cross section based on such a gauge- dependent amplitude (defined 
in the 't Hooft-Feynman gauge), known as "CC03 cross section", is a very useful quan- 
tity that is very close to the full 4/ calculation if both W bosons are close to resonance. 
The CC03 amplitude can be rendered gauge invariant upon deforming the momenta of 
the four outgoing fermions in such a way that the intermediate W-boson states become 
on shell, because the residues of the W resonances are gauge-invariant quantities. This 
"on-shell projection" , which is part of the DPA, involves some freedom, and different ver- 
sions, which have been described in Refs. [15,17], differ by contributions of relative order 

e>(r w /M w ). 
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We want to perform the exercise to study the usefulness of a possible "CC03" 3 off- 
shell cross section for 77 — > WW — > 4/. To this end, we define the amplitude for the 
off-shell W-pair signal by evaluating the three W-pair diagrams in the non-linear gauge 
with polarization vectors Si{ki) for the incoming photons, which obey the gauge conditions 



£i(fci) • k 2 = e 2 {k 2 ) ■ h — 0. 



(2.28) 



In terms of WvdW spinors, this means that the gauge spinors g\ and g 2 of the photons 
are identified with the spinors of the momenta k 2 and A7, respectively. With this choice 
the auxiliary functions for the matrix elements (2.3) read 



A ++,ww( k i,Pj,Qj) = Pw(pi +P2)Pw(P3+P4 



(P2P4)* 

(hk 2 y 



X 



Pwipi +P2- h) 

(pm) 



+ 



{hk 2 y 



(k 2Pl )(k 2 p 3 )(k 2 [P 1 + p 2 ]h) + (k lPl )(k lP3 )(k 1 [P 3 + P 4 ]k 2 ) 

(UP! + P 2 ]k 1 )(k 1 [P 3 + Pt\k 2 ) - 2(h • k 2 ) (k m ) (k 2 p 3 ) 



- 2 (PiPs) (^1^2) J + (h <-> fc 2 ) , 

^+-,Ww(^Pj 5 Qj) = Av(pi +P2)-Pw(P3 +Pi) 



X ^W(Pl + P2 - fcl) + PwiPl + P 2- k 2 ) 



r w =o 



x |(fc 2 [Pi + P 2 ]A;i 

_ (P2P4)* (pm) 

2(hk 2 ) 



( hp 2 ) * ( k 2 p 4 ) * (pip 3 ) (p 2 p 4 )*(k 1 p 1 )(k 1 p 3 ) 



(k\k 2 



(k\k 2 



(k 2 [Pi + P 2 ]h) 2 + (k 2 p 2 )*(k m )*(k lPl )(hp 3 ) 



(2.29) 



Note that ^a^a^ww d° no ^ coincide with the parts of the functions ^AiA2,w °f ^q. (2.10) 
that are proportional to Pw(pi + p 2 )-P\v(P3 + Pa) because the derivation of Eq. (2.10) 
involves rearrangements of various singly-resonant contributions. We point out that the 
definition (2.29) is neither independent of the gauge fixing used to define gauge-boson 
propagators nor of the gauge of the external photons. The definition is gauge invariant 
after the outgoing fermion momenta Pi are on-shell projected as described above, while 
leaving the resonant propagators Pwipi + P2)Pw(P3 + Pa) untouched. This defines the 
lowest-order amplitude in DPA. Finally, we stress that the t- and w-channel W propagators 
in Eq. (2.29) do not receive a finite W width; otherwise the gauge invariance of the DPA 
would be spoiled. 



3 Inclusion of anomalous gauge-boson couplings 

In this section we introduce the most important anomalous gauge-boson couplings 
accessible by the process 77 — > 4/ and give explicit analytical results for the corresponding 
helicity amplitudes. 

3 The name also fits to the 77 case where three W-pair diagrams exist in unitary or non-linear gauge. 
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3.1 The effective Lagrangians 



First we consider anomalous triple gauge-boson couplings (ATGC) in the charged- 
current sector, i.e., anomalous 7WW and the related 77WW couplings. Instead of using 
rather general parametrizations of non-standard couplings [29], we follow the approach 
already used at LEP2 to reduce the number of free parameters by requiring that all 
symmetries of the SM are respected. From the resulting operators we only keep those 
that appear in the lowest-order cross section of 77 — > 4/. Specifically, we start from the 
gauge-invariant CP-conserving effective Lagrangian with dimension-6 operators [30] 



aw 



(3.1) 



where $ is the Higgs doublet field and 



= d^B v - d u B>*, 
W"" = (WF, Wf, Wf) = d"W - d"W» + g 2 W x W 



(3.2) 



are the field strengths of the U(l) and SU(2) gauge fields, respectively. The Pauli matrices 
are combined into the vector a = (a±, a 2 , 03), and the parameters g±, g 2 denote the gauge 
couplings. 4 Inserting the vacuum expectation value of the Higgs field $, we can relate 
the coefficients 0^0, ®w<f>, an d «w to the coefficients of the Lagrangian considered in the 
LEP2 analysis [30], 



Ag\ 



«W0 



■(Ak z - Agf) = a W(j} + a B( p, A 7 = A z 



(3.3) 



In contrast to the pure anomalous 7WW coupling [29], the SU(2)xU(l) symmetry of the 
effective Lagrangian (3.1) induces additional anomalous 77 WW and 7W0 couplings. The 
corresponding Feynman rules are 

ir 7W+w ^ k _j = _ ie | A/c^fcoi/^v - ko P 9i^) 

k +ll k^ v k 0p - k- p k +p k 0v + g^k-^k+ko) - k +p (k_k )) 



+ g w (k v(k+ k -) ~ k_ v (k + k )) + g^ik+pik-ko) - k 0p (k + k_)) 



ir 



77W+W- 



\ivpa 



(fa, k 2 , k+, k-) 



-le 



2 A 7 



guvg P a(ki + k 2 ) 2 + g PP g ua (k 2 k + + fafa) 



+ 9u P g^a{k 1 k + + k 2 k-) + g pv (fa + k 2 ) p k +a + (fa + k 2 ) a k„ p 
+ 9pa (k+ + k_) p fa v + (k + + k_) v k 2p + g pp (fa — k 2 ) a k +v — fa u k +a — fa a k_ v 



4 In order to be compatible with the conventions of Ref. [36] used for the SM amplitudes above, we 
had to change the sign of the SU(2) coupling g 2 w.r.t. Ref. [30]. 
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+ 9/j.a {ki — k2) p k- v — k\ v k- p — k\ p k +v + g vp {k,2 — ki) a k +pj — k2^k +a — k2 a k-p 



+ 9v 



(k,2 — ki) p k_^ — k2f 1 k_ p — k 2p k + p 



LS.K f 1 

ir^(/c , fc w , hf) = -iej^-^k^kojg^ - fc^fc 0>v |, (3.4) 

where all fields and momenta are considered incoming. Note that the neglect of the 
contribution to the quartic coupling 77WW, which is proportional to A 7 , would lead to 
a violation of electromagnetic gauge invariance in predictions for 77 — > WW(^ 4/). In 
contrast, neglecting the 7W0 coupling, which is proportional to Ak 7 , would not spoil the 
electromagnetic gauge invariance of the predictions. 

Next we consider anomalous triple gauge-boson couplings involving only the neutral 
gauge bosons 7 and Z. Assuming Lorentz invariance and electromagnetic gauge invariance, 
the most general effective dimension-6 Lagrangian for 77Z, 7ZZ, and ZZZ couplings can 
be written as [43] 5 

£nc gc = ( v) - Hidu.znw + [/?( v) - fi(d,zn\Kz p 

+ [hJ(dpFn - h\{d,Z^)\F vp Z p + [hJ(dpFn - hl(dpZn]F up Z^ (3.5) 

with the abelian field-strength tensors 

F pv = d"A u - d v A», Z^ = d p Z v - d v Z p (3.6) 
and the dual field-strength tensors (e 0123 = +1) 

\ iivpc rp y uv _ 1 

An operator inducing a 777 coupling does not appear in Eq. (3.5) since it violates elec- 
tromagnetic gauge invariance. 

Apart from the 77WW coupling which is induced by symmetries in the Lagrangian 
(3.1), we also include genuine anomalous quartic gauge-boson couplings (AQGC) in our 
analysis, whose lowest dimension is 6. In Refs. [31,32] all genuine dimension-6 AQGC 
that involve photons and that are allowed by electromagnetic gauge invariance and cus- 
todial SU(2) C have been classified; more general AQGC have been discussed in Ref. [44]. 
Following Ref. [32] we use the effective Lagrangian 

2 

C i?vv = - F ^ F ^W Q W Q + a c F^F^W a + ~a F^W Q W a } (3.8) 

with the definition 

W, = (Wl,Wl,Wl) = + W-)„ -±=(W + - W-)„ i-Z„) . (3.9) 

The scale of new physics, A, is introduced in Eq. (3.8) to render the coupling coefficients 
do, a c , ao dimensionless. The effective Lagrangian 

C^vv contains 77WW and 77ZZ cou- 
plings, whose Feynman rules can be found in Ref. [32] . The other coupling structures C n 
and C n considered in Ref. [32] induce 7ZWW couplings that are not relevant for 77 — > 4/. 

5 Note that our conventions differ from those of Ref. [43] by a minus sign in the Z-boson field. 



F^ = -e^ pa F p(T , Z*" = -e pvp °Z pa . (3.7) 
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(a) (b) (c) 




Figure 2: Representative diagrams with anomalous 7WW and 77WW couplings (black 
blobs) contributing to CC processes 77 — > 4/. 

3.2 Amplitudes with triple gauge-boson couplings 

Before we write down the helicity amplitudes including ATGC explicitly, we discuss the 
impact of these couplings w.r.t. the SM cross section. The diagrams containing ATGC and 
the corresponding quartic couplings in CC diagrams are shown in Figure 2. We quantify 
the size of the anomalous contributions in terms of powers of anomalous coupling factors 
(generically denoted by a 3 ) or suppression factors r w /M w . Considering the SM process 
77 — > WW — > 4/ as the leading contribution, i.e., regarding anomalous-coupling effects 
as small, we get non-standard contributions to CC and CC/NC cross sections from CC 
ATGC of the following orders: 

• 0(a 3 ): 

The matrix elements of diagrams (a) and (b) in Figure 2 involve one power of 03. 
Both diagrams are not suppressed by T w /M w since they are doubly resonant. 

• 0(a 3 T w /M w ): 

The diagram (c) of Figure 2 has one power of a 3 and one resonant W-boson propa- 
gator, i.e., it is only singly resonant. Thus, it is of C(a 3 r w /M w ). 

• 0{4): 

The diagrams (d) and (e) of Figure 2 involve two anomalous couplings a 3 and are 
doubly resonant. Therefore, they are of C(a|). Note that the squares of the diagrams 
(a) and (b), as well as their products with one another, are of the same order as the 
interference of diagrams (d) and (e) with the SM amplitude. 

There are no diagrams containing CC ATGC for NC processes. 
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Next we consider the impact of NC ATGC, as defined in the effective Lagrangian (3.5). 
The by far largest SM cross sections of the process class 77 — > 4/ belong to diagrams 
with two resonant W bosons in CC and CC/NC reactions. Thus, the largest effect of 
NC ATGC could be expected from an interference of "anomalous diagrams" with the SM 
amplitude for CC or CC/NC processes. The only candidate of this kind is a diagram 
where an off-shell s-channel Z boson is produced by an anomalous 77Z coupling that 
subsequently produces a W-boson pair. However, the effective 77Z coupling of Eq. (3.5) 
vanishes for two on-shell photons, so that this diagram does not contribute. No other CC 
diagram exists that includes a NC ATGC. 

We now turn to the effects of NC ATGC in NC amplitudes, i.e., in diagrams without 
W bosons. The corresponding SM amplitudes involve at most a single resonance of the 
Z boson, which leads already to a suppression of NC cross sections w.r.t. CC cross sections 
by a factor (Tz/Mz) 2 - This suppression is clearly visible in the numerical results presented 
in Section 6.2. f below. Diagrams with one NC ATGC also possess at most one resonant 
Z boson and, therefore, show a suppression by a factor a^iYz/ Mz) 2 w.r.t. the CC signal 
diagrams. This suppression is not changed by interferences with doubly-resonant CC 
diagrams in CC/NC processes because the Z- and W-boson resonances are located at 
different regions in phase space and do not enhance each other. Diagrams with two NC 
ATGC can involve two Z-boson resonances resulting in a suppression of 0(alTz/Mz), 
which is also small compared to the CC case owing to the squared ATGC. In summary, 
we conclude that the sensitivity of the processes 77 — > 4/ to NC ATGC is much smaller 
than to CC ATGC. Therefore, we restrict our investigation on ATGC to CC couplings in 
the following. 

As explained above, the diagrams of Figure 2 induce contributions to the amplitude 
that are either linear or quadratic in the CC ATGC. We give the explicit contributions 
to the helicity amplitudes in a way similar to the SM case (2.3), 



Mx\xl°ccATGc(kh Pj, Qj) — e4 ^o- 1 ,-5 -2,+^o-3,-^cr4,+ 5 , w / 1 j 2 fi , w / 3 j 4 ^3Ax\xl(ki,pj, Qj) (3.10) 



amplitudes for the various process types proceeds as described in Sections 2.3.2 and 2.3.3. 
The terms in <5 3 A^^ that are quadratic and linear in ATGC explicitly read 





quad 



Pw(pi +P2)Pw{P3 +Pi)Pw(pi +P2- ki){k 1 p 1 ){k 2 pz) 




+ (PiP2}*(P3P4)*((kiP3}{k 2 Pi} - {k l k 2 ){pm)) 



1 



(P1P2) * {PzPa) *{k\Pi ) {k 2 p 3 ) 



+ (kik 2 )({p 3 p 4 )* (hp^* (kipz) - (pip 2 )*{k 2 P4)*(k 2 Pi)) 



^ttHpiP2)*(P3P4)*^(pi +P2 - ki) 2 {k 1 p 3 ){k 2 p 1 ) - {hk 2 ) {p^) 



w 



+ {h<r+k 2 ), 



Mil 



quad 



Pw(pi + P2)Pw(P3 + P4)Pw(pi +P2- ki){k 2 p A )*{kipi) 



15 



X -Ak 2 



w 



+ A/c 



7 M 2 



-2(pi +p 2 - h) 2 (hP2)* (km) 



+ ( P2 [K 2 - KMfaiPi + p 2 ]h) - ( PlP2 y{k 2P4 y{ P3P4 ){k lPl ) 



+ 



A? 



{piP2)*{pm) 



~2^ Pl +P2 ~ kj 2 (k 2 p4)* (hpi) 



+ <pt[K 2 -P 3 ]k 1 )(k 2 [K 1 -P 2 ] Pl ) 

+ ({pi, Qi,P2, Q2} ^{ps, Qs; pa, Qa}) , 



6*A; 



X 



llin 



2-Pw(Pi +P2)Pw(P3 +pd(k 2P i)(k 2 p 3 ) 



{P2PaY 



A/t 7 (j9 2 p 4 ) H 



A, 



(P3P4)*(PlP2)*{PlP3) 



{k lP2 )*(k lP4 )* 

+ |2(Q 4 - <2 3 ) [-Qi + (Qi - Q2)2(k lPl )P w ( Pl + pa)] P W (P3 + Pa) 

(k2P 3 }(P2[Pl-K 1 }k 2 ) 



X Pw(P3 +P4 - k 2 )- 



(k lPl )*(k lP2 y 



x 



A 



S 3 A: 



x 



A«: 7 (p 2 p 4 )* + — y Y (pm)*{p2[P* - K 2 } P3 ) 
+ ({pi,Qi;p2,Q2}<->{p3,Q3;p4,Q4}) f + (fci 



-A/t 7 (p 2 p 4 )* + -^-(( P i P2 y{k 2P 4y{k 2Pl ) - {pzPa)* {k 2 p 2 )* {k 2 p 3 ) 
+ (P1P2}* (PsPa)* {PiPs}) 

2(Qi - Q 2 )Pw(pi + P2) [[Q 4 + (Q 3 - Q 4 )2(A; 2 p 4 )Pw(P3 + p 4 )] 



x 



+ 



A. 



{PiP2)*(piP3)(P3+Pa - k 2 ) 



+ ({pi,<5i;p2,<52>^{p3,<53;p4,<54})| 

+ (c.c. and {p ls Qi;p 3 , (? 3 ; &i} ^{P2, Q^Pa-, Qa, k 2 }), 



(3.11) 
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Figure 3: Diagram with AQGC (black blob) contributing to 77 — > 4/. 



where "c.c. and {. . .} «-►{. . .}" indicates that the complex conjugate of the preceding 
expression has to be added after some substitutions. The auxiliary functions for the 
remaining polarizations follow from the relations (2.11) and (2.12). 

In order to check our results, we have implemented the ATGC of the effective La- 
grangian (3.1) into the program Madgraph [34] and compared our amplitudes with the 
Madgraph results for various phase-space points. We found perfect numerical agree- 
ment. 



3.3 Amplitudes with genuine quartic gauge-boson couplings 

Figure 3 shows the only diagram with an AQGC (generically denoted by 04) that 
contributes to 77 — > 4f. For CC processes the "anomalous diagram" contributes in 
C(a 4 ) to the cross section, because it is (as the SM contribution) doubly resonant. For 
NC processes, the diagram involves one power of 04 and two Z-boson resonances and 
interferes with the singly-resonant SM amplitude. In this case, the contribution to the 
corresponding cross section is suppressed by aJPz/Mz w.r.t. CC cross sections, i.e., the 
suppression factor involves one factor in the anomalous coupling or in Tz/Mz less than we 
counted for NC ATGC. In the following we take both CC and NC AQGC into account. 

The AQGC contributions to the amplitudes read 

e 4 

Mllll^v = —5 ai ^ a2 5 a3 ^ a4 g^ vv gyy ih gy 3 hh P v (p 1 +p 2 )Pv(p3+P4) 

x8 A AllZ(k l ,k 2 ,p 1 ,p 2 ,p 3 ,p A ) (3.12) 



with 



and 



1 

3rrww — 1, SVyZZ — — (3.13) 



hA ++ (ki,k 2 ,Pi,P2,P3,P4) = (4a - 4ia + a c ){p 2 p A )*{k 1 k 2 ) 2 {p 1 p 3 ), 
5 4 A~Z(k 1 ,k 2 ,p 1 ,p 2 ,p 3 ,p 4 ) = -2a c (k 2 p 2 )*(k 2 piY(kxpi)(kip z ). (3.14) 

The remaining auxiliary functions S^A^^ can be obtained via the substitutions 

^A a x 1 1 ' x + 2 (k 1 ,k 2 ,p 1 ,p 2 ,p 3 ,p 4 ) = hAH'x^k^k^p^p^p^ps), 
hA^;l 3 2 (k u k 2 ,p u p 2 ,p 3 ,p 4 ) = ^^(fcij^.pa.pi.ps,^), 
hAl 1 l a x s 2 (k 1 ,k 2 ,p 1 ,p 2 ,p 3 ,p4) = (hAZ' J x 1 1 ]Zx i 2 (k 1 ,k 2 ,p 1 ,p 2 ,p 3 ,p 4 )^ . (3.15) 
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{W,Z} 



H 




{W,Z} 



Figure 4: Diagram with effective 77H coupling (black blob). 



The generic amplitude M-\^\^Zvv is coherently added to the SM amplitude M." 10 1 



of Eq. (2.3) for V = W, Z, respectively. The colour summation of the squared amplitudes 
for the various process types proceeds as in the SM case. 

Again we have checked the amplitudes against results obtained with Madgraph, as 
explained at the end of the previous section. 



4 Effective 77H coupling and Higgs resonance 

In order to incorporate a possible Higgs resonance in 77 — > H — > VV — > 4/ with 
V = W, Z, as depicted in Figure 4, we consider an effective coupling of the Higgs boson 
to two photons. In the SM this coupling is mediated via fermion (mainly top-quark) and 
W-boson loops. We define the effective Lagrangian for the 77H vertex [45] by 



77H 



Pf" TP 

- 1 - 1 fll/ 1 



H 

V 



(4.1) 



where v = 2Mws w /e is the vacuum expectation value of the Higgs field H. Up to 
normalization, £ 77 h is the lowest-dimensional, CP-conserving, electromagnetically gauge- 
invariant operator for two photons and the scalar field H. The corresponding Feynman 
rule reads 

ir^ H (A; 1) k 2 , fen) = ^ [g^ihh) - k^k u 



<>MJ ' 



(4.2) 



where ki, k 2 are the incoming photon momenta. Comparing this Feynman rule to the 
loop-induced SM vertex with the external fields on shell, which has, e.g., been given in 
Refs. [20,45], we obtain 



w 



SM 



7T 



Ml 



+ 1 + 



6M 



w 



2 E N fQ 

f 



Ml 
2 m) 



(2M^-Ml)C (M R ,M, 



w 



2 + {4m}-Ml)C {M H ,m f ) 



(4.3) 



where the colour factor iVS in the sum over all fermions / is equal to 3 for quarks and 1 
for leptons. The scalar 3-point integral Cq is given by 



C (M H ,m) 



2Ml 



In 2 



'Pr, 



1 



A, 



1 



4m 2 
Ml 



+ i0. 



(4.4) 
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The complete matrix elements for the diagrams with a Higgs resonance (as shown in 
Figure 4) can then be written as 

e 4 

^xlxl^HVV = ~7TY K,-aJa 3 ,-a 4 £ 77 H 9^VV 0y/ 1 / 2 0y/ 3 / 4 i V(Pl + P2)Pv(P3 + Pi) 

x M^P u (k 1 + k 2 )5nAll%(k 1 ,k 2 ,p 1 ,p 2 ,p 3 ,p 4 ) (4.5) 
with g-y 7 vv defined in Eq. (3.13) and 

5nA--(k 1 ,k 2 , Pl ,p 2 , P3 ,p 4 ) = ^M_(p 2 p 4 y(p lP3 ), ShA^ = 0. (4.6) 

The other expressions for S^A^J^ follow in the same way as described in Eq. (3.15) for 
S^A^^. The width in the Higgs-boson propagator P H is introduced in the same way as 
in Section 2.4 for the gauge bosons. 



5 Monte Carlo integration 

The squared matrix element is integrated over the phase space following the strategy 
described in Refs. [16,41,46]. We apply the multi-channel Monte Carlo technique [25] 
where the integrand is flattened by choosing appropriate mappings of the pseudo-random 
numbers into the momenta of the outgoing particles. This is necessary because the inte- 
grand shows a very complex peaking structure in eight (77 — > 4/) or eleven (77 — > 4/7) 
dimensions. This complicated structure of the integrand (induced by various diagram 
types) deteriorates the statistical error of the numerical integration and can lead to nu- 
merically unstable results. For each diagram we introduce an own appropriate mapping, 
so-called "channels" , so that more events are generated in regions where the squared ma- 
trix element of the diagram is large. The multi-channel approach combines these channels 
in such a way that the whole integrand is widely smoothed everywhere. For each event 
a single channel is randomly chosen, and the phase-space configuration is determined ac- 
cording to the mapping of this channel. The probability to choose a given channel, called 
a-priori weight, is optimized according to Ref. [26] to minimize the statistical error as 
much as possible. 

The convolution over the photon spectrum is given by 

do- = J d^i^ dx 2 f 1 (x 1 )f 1 (x 2 )da JJ (x 1 P 1 ,x 2 P 2 ), (5.1) 

where da 77 is the differential 77 cross section. The function / 7 (xj) denotes the probability 
density for getting a photon with momentum ki = XiPi, and Pj is the electron momentum 
before Compton backscattering. In order to reduce the statistical error of this integration 
we use a simple way of stratified sampling. The integration region for Xi of each photon 
spectrum is divided into a fixed number of bins. We choose bin i with a probability «j 
and divide the corresponding weight by c^. In this way the integration remains formally 
unchanged if we normalize J2i a i — 1- The parameters ctj can be used to improve the 
convergence of the numerical integration. By choosing the ctj proportional to the cross 
section of the corresponding bin i, more events are sampled in regions where the photon 
spectrum is large. Care has to be taken that the ctj do not become too small because this 
might lead to rare events with very large weights that render the error estimate unreliable. 
This optimization typically reduces the Monte Carlo integration error by a factor 2-5. 
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6 Numerical results 



6.1 Input parameters 

We use the following set of input parameters [47]: 



M w = 


80.423 GeV, 


r w — 


2.118 GeV, 


M z = 


91.1876 GeV, 


r z = 


2.4952 GeV, 


M H = 


170 GeV, 


r H = 


0.3834 GeV, 


a(0) = 


1/137.03599976, 


a s = 


1.1172, 




1.16639 x 10- 5 GeV~ 2 , 







(6.1) 

where the Higgs mass is chosen well above the W-pair threshold so that intermediate 
Higgs bosons decay rapidly into W pairs; the corresponding decay width Th has been 
obtained with the program HDECAY [48]. 
Furthermore, we apply the separation cuts 

£ 7 >10GeV, 0(7, beam) > 5°, 9(1, 7) > 5°, %,7)>5°, 

Ei> 10 GeV, #(/,beam)> 5°, 9(1,1') > 5°, 9(l,q)>5°, 

E q > 10 GeV, 9(q, beam) > 5°, m(q, q') > 10 GeV, (6.2) 

where q and / denote quarks and charged leptons, respectively, and m(q, q') is the invariant 
mass of an outgoing quark pair. The energies Ex and angles 9(X, Y) are defined in the 
laboratory frame. Using these cuts all infrared, i.e., soft or collinear, singularities are 
removed from the phase space. 

In order to account for leading universal corrections, we use two different values for 
the coupling constant a = e 2 /(An). Since on-shell photons couple to charged particles 
with the coupling constant a(0) (effective electromagnetic coupling at zero-momentum 
transfer), we take this coupling for each external photon in the processes 77 — * 4/ and 
77 — > 4/7. For CC reactions, the remaining couplings correspond to W// vertices. For 
these vertices a large part of the electroweak radiative corrections [23] (the running of the 
electromagnetic coupling and the universal corrections related to the p parameter) are 
absorbed into an effective electromagnetic coupling a^^ which is derived from the Fermi 
constant by 

7T 

Therefore, in the following numerical studies, we replace a 4 by a(0) 2 o;^ for the processes 
77 — * 4/ and a 5 by a(0) 3 aQ^ for 77 — > 4/7. 

For the evaluation of the photon spectrum we use the program CompAZ [27] with 
the polarization of the laser beams —1 (i.e. photon helicity — l) 6 and the polarization of 
the electron beams +0.85. This choice for the relative signs in the polarizations yields 
a sharper peak at the upper end of the photon spectrum. Results for monochromatic 
photon beams are always shown for unpolarized photons. 

6 Internally in CompAZ the polarization of the laser light is defined as the negative of the photon 
helicity. 
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If not stated otherwise, the results are obtained in the fixed-width scheme. 

The numerical integration is carried out using 10 7 events. The runtime of our Monte 
Carlo program on a PC with 2 GHz varies from 30 minutes to 6 hours depending on the 
considered process. 

6.2 Results for integrated cross sections 

6.2.1 Survey of cross sections 

In order to illustrate the reliability of our Monte Carlo generator we compare our 
results on cross sections for a representative set of the processes 77 — > 4/ and 77 — ■> 4/7 
with the results obtained with the Monte Carlo program Whizard (version 1.28) [33] 
which uses the matrix-element generator Madgraph [34] 7 . In Tables 2 and 3 we list the 
results for the 17 different final states defined in Table 1. The numbers in parentheses 
correspond to the Monte Carlo error. For the final states that can be produced via inter- 
mediate gluons we compute the cross section both with and without gluon-exchange con- 
tributions. Since the version of Madgraph implemented in Whizard is not able to deal 
with interferences of electroweak and QCD diagrams, we give only the pure electroweak 
Whizard/Madgraph results for these processes. Furthermore, we list the correspond- 
ing cross sections with and without convolution over the photon beam spectrum. For this 
study, we have implemented the program CompAZ into Whizard. 

As explained in Section 2.2, the cross sections for the CP-equivalent final states 
e^u c ud( , ~f) and u e e + du('j) are not identical if the convolution over the photon beam spec- 
trum is carried out. Therefore, we give results for both final states. In all other cases, the 
cross sections for a given final state and for the CP-conjugated one coincide. 

CC and CC/NC processes possess the largest cross sections because of the dominance 
of W-pair production. The convolution over the photon spectrum reduces these cross 
sections significantly since low-energy photons cannot produce on-shell W pairs. NC 
processes are affected less, and in some cases, such as 77 — > e + e~/i + /i~, the cross section 
is only slightly reduced. Owing to the colour factors of the quarks, hadronic and semi- 
leptonic cross sections differ by roughly a factor 3, hadronic and leptonic cross sections 
by roughly a factor 3 2 = 9. For CC processes 77 — > 4/ we obtain a rough estimate 
of the cross sections by multiplying the cross section of 77 — > WW with the branching 
ratios of the W bosons into leptons or quarks depending on the final state. Note that this 
estimate, which is only good within 10—20%, does not take into account contributions 
from background diagrams, width effects, and cuts on final-state fermions. The difference 
of cross sections for CC processes and the corresponding processes of mixed type reflects 
the size of the background contributions induced by NC diagrams. 

The results of Whizard, which are also generated with 10 7 events, and of our program 
typically agree within 1-2 standard errors. The size of the statistical errors obtained with 
Whizard and our program is comparable. The runtime of Whizard is usually somewhat 
bigger than the one of our program. Depending on the process class, the speed of our 
program is 1—7 times higher, where the largest difference occurs for NC processes. 

7 For a tuned comparison we rescaled the Whizard/Madgraph results by a factor a(0) 2 a^/a 4 for 
77 — * 4/ and a(0) 3 «g /a 5 for 77 — > 4/7. 
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77 -> 


present work 


Whizard /Madgraph 


(T4/ [ fb] 


(74/ [ fb] (conv) 


(7 4 / [fb] 


(74/ [fb] (conv) 




826.47(21) 


190.87(10) 


826.39(26) 


191.05(16) 




1.75460(62) 


0.90525(61) 


1.75518(78) 


0.9050(11) 




19.400(33) 


19.129(61) 


19.342(21) 


19.188(48) 


e~e + e~e + 


9.469(17) 


9.357(32) 


9.453(11) 


9.383(25) 


e~e + z/ e z/ e 


828.34(21) 


191.72(10) 


828.29(26) 


191.55(17) 


e~z/ e ud 


2351.11(68) 


565.05(33) 


2351.79(84) 


565.07(51) 


z/ e e + du 


2350.84(68) 


558.39(32) 


2353.21(84) 


558.41(50) 


z/ c z/ e uu 


1.19761(50) 


0.61256(50) 


1.19684(57) 


0.61083(71) 


z/ e z/ e dd 


0.095981(44) 


0.049092(45) 


0.096011(48) 


0.049118(57) 


e~e + uu 


14.036(21) 


10.597(26) 


14.016(15) 


10.574(21) 


e~e + dd 


4.7406(29) 


2.6614(32) 


4.7377(28) 


2.6651(38) 


udsc 


6659.6(2.1) 


1603.8(1.0) 


6663.5(2.7) 


1605.0(1.5) 


uucc 


10.469(14) 


6.111(12) 


10.4531(88) 


6.113(10) 


with QCD 


1543.6(2.9) 


1071.3(2.9) 


— 


— 


uuss 


3.3282(21) 


1.6569(18) 


3.3310(20) 


1.6595(23) 


with QCD 


412.97(75) 


288.79(72) 


— 


— 


ddss 


0.49807(29) 


0.23232(24) 


0.49804(30) 


0.23252(32) 


with QCD 


96.34(18) 


66.80(18) 






uuuu 


5.1846(69) 


3.0298(57) 


5.1900(45) 


3.0419(53) 


with QCD 


772.6(1.5) 


538.9(1.4) 






dddd 


0.24683(15) 


0.11581(12) 


0.24665(17) 


0.11579(17) 


with QCD 


48.252(96) 


33.685(88) 






uudd 


6663.5(2.3) 


1606.1(1.1) 


6664.8(2.8) 


1604.6(1.6) 


with QCD 


7075.8(3.7) 


1896.4(2.9) 







Table 2: Total cross sections for 77 — > 4/ at y/s = 500 GeV for various final states with 
and without convolution over the photon spectrum. 
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77 -»• 


present work 


Whizard /Madgraph 


(74/7 [fb] 


<74/ 7 [ fb] (conv) 


<7 4/7 [fb] 


04/ 7 [fb](conv) 


e-z/ e z/ M /i+7 


39.234(44) 


6.188(11) 


39.218(29) 


6.2040(87) 




0.10157(10) 


0.028612(40) 


0.101556(88) 


0.028548(52) 


e~e + /i~/i + 7 


1.0567(35) 


0.5083(28) 


1.0547(20) 


0.5091(29) 


e~e + e~e + 7 


0.5085(18) 


0.2433(13) 


0.5091(10) 


0.2461(12) 


e"e + i/ e P e 7 


39.301(46) 


6.213(11) 


39.332(30) 


6.2069(89) 


e~z/ e ud7 


96.61(13) 


14.216(27) 


96.575(75) 


14.159(21) 


z/ c e + du7 


96.60(13) 


15.459(30) 


96.520(76) 


15.429(22) 


z/ e z/ e uu7 


0.030818(35) 


0.008640(14) 


0.030756(28) 


0.008609(16) 


z/ c z/ c dd7 


0.00061753(75) 


0.00017313(31) 


0.00061731(56) 


0.00017358(34) 


e~e + uu7 


0.6446(17) 


0.25463(99) 


0.6477(10) 


0.2579(10) 


e~e + dd7 


0.26653(36) 


0.08137(17) 


0.26689(28) 


0.08166(21) 


udsc7 


229.86(36) 


32.621(81) 


229.52(19) 


32.531(49) 


uucc7 


0.30556(69) 


0.10718(34) 


0.30563(47) 


0.10836(43) 


with QCD 


34.73(14) 


13.801(77) 


— 


— 


uuss7 


0.08791(13) 


0.026278(59) 


0.087935(98) 


0.026271(65) 


with QCD 


6.362(23) 


2.493(13) 


— 


— 


ddss7 


0.0046253(71) 


0.0014842(37) 


0.0046191(52) 


0.0014832(36) 


with QCD 


0.5427(22) 


0.2165(11) 






uuuu7 


0.15081(33) 


0.05301(16) 


0.15082(21) 


0.05332(16) 


with QCD 


17.377(71) 


6.964(35) 






dddd7 


0.0022893(37) 


0.0007421(21) 


0.0022878(25) 


0.0007398(18) 


with QCD 


0.2716(11) 


0.10863(53) 






uudd7 


229.86(40) 


32.85(15) 


229.65(19) 


32.518(51) 


with QCD 


236.31(42) 


35.14(11) 







Table 3: Total cross sections for 77 — > 4/7 at y^s = 500 GeV for various final states with 
and without convolution over the photon spectrum. 
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Figure 5: Integrated cross sections of the processes 77 — ► e z/ e ud(7) with and without 
convolution over the photon spectrum as a function of the CM energy y's. 



6.2.2 Energy dependence of integrated cross sections 

In Figure 5 we show the cross sections for the processes 77 — > e~u e nd( , j) as a function 
of the centre-of-mass (CM) energy y/s with and without convolution over the photon 
spectrum. Here and in the following, with convolution over the photon spectrum yfs 
stands for the CM energy t/s^ of the incoming electron beams, without convolution it 
is the CM energy ^/s 77 of the incoming photons. In the case without photon spectrum, 
the rise of the cross section is clearly visible at the W-pair threshold, y/s^ > 160 GeV. 

For 77 — > e~z/ e ud the cross section increases roughly proportional to /3 = ^/l — 4M^/s 77 
above the threshold, as expected from the two-particle phase space of the W pairs. For 
77 — > e~u e VLd'j the rise of the cross section is not as steep because of the higher-dimensional 
WW7 phase space. The convolution over the photon spectrum reduces the available 
energy for W-pair production and shifts the onset of the cross section to higher CM 
energies. 

The cross sections for 77 — > 4/ as well as 77 — > Af'j decrease at high energies, even 
though the total cross section of the 77 — > WW process approaches a constant in the 
high-energy limit if no cuts are imposed, i.e., if the W bosons are allowed to go in the 
beam directions. At high energies, however, forward and backward scattering of W bosons 
is restricted due to the cuts applied to the outgoing fermions, because the decay fermions 
mainly follow the direction of the decaying W boson. 

6.2.3 Contributions from CC, NC, and gluon- exchange diagrams 

In Figure 6 we show the impact of CC, NC, and gluon-exchange diagrams on the 
CC/NC processes 77 — > uudd and 77 — > uudd7. We do not include the photon spectrum 
in this analysis. Above the W-pair threshold, Js^ > 160 GeV, the CC diagrams are 
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Figure 6: Different contributions to the integrated cross sections for the processes 77 
uudd (7) as a function of the CM energy without photon spectrum. 



clearly dominating, while the contributions from gluon-exchange diagrams are one or 
two orders of magnitude smaller. The impact of the gluon-exchange diagrams strongly 
depends on the choice of the invariant-mass cut between two quarks, and gluon-exchange 
diagrams are more important if the invariant-mass cut is small. The contributions from 
pure NC diagrams are totally negligible as long as W-pair production is possible. 

6.2.4 W-pair signal diagrams and double-pole approximation 

In Figure 7 the cross sections of the W-pair signal diagrams and the DPA for 77 — > 
WW — > 4/ (see Section 2.5 for definitions) are compared with the complete lowest-order 
cross section for several processes. The plots on the l.h.s. show the cross sections for 
various final states calculated from the full set of (electroweak) diagrams, from the signal 
diagrams only, and in DPA separately for hadronic, semi-leptonic, and leptonic final 
states, while the plots on the r.h.s. show the relative deviation from the corresponding 
DPA. We do not include the convolution over the photon spectrum and gluon-exchange 
diagrams in this analysis so that effects of the approximation are clearly visible. For 
energies not too close to the W-pair threshold, the DPA agrees with the full lowest-order 
cross section within 1-3%, which is of the expected order of r w /M w . Near threshold, 
i.e. for y/s^ — 2Mw = 0(Tyr), the reliability of the DPA breaks down, since background 

diagrams become more and more important and small scales 7, such as yj s 77 — 4M^, can 
increase the naive error estimate from r^/Mw to Tw/7- The cross sections of the W-pair 
signal diagrams, however, shows large deviations from the full 77 — > 4/ cross sections for 
the whole energy range, in particular, at high energies. As explained in Section 2.5, the 
W-pair signal diagrams are not gauge invariant, and thus the reliability and usefulness of 
the resulting predictions should be investigated carefully. The results of Figure 7 clearly 
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Figure 7: Cross sections of various processes including all diagrams, only W-pair signal 
diagrams, and in DPA as a function of the CM energy (l.h.s.), and the corresponding 
relative deviations from the DPA (r.h.s.); photon spectrum and gluon-exchange diagrams 
are not included. 
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<r(77 -»• e ^c^/U + ) 


[ GeV] 


500 


800 


1000 


2000 


10000 


fixed width 


826.40(21) 


788.35(21) 


746.94(21) 


500.70(20) 


31.745(68) 


step width 


827.45(22) 


789.34(21) 


748.17(23) 


501.41(21) 


31.746(68) 


running width 


827.43(23) 


789.29(21) 


748.11(23) 


501.32(21) 


31.715(68) 


complex mass 


826.23(21) 


788.18(21) 


746.78(21) 


500.59(20) 


31.738(68) 




(7(77 -> e z/ c ^ + 7) 


v /s^ [ GeV] 


500 


800 


1000 


2000 


10000 


fixed width 


39.230(45) 


47.740(73) 


49.781(91) 


43.98(18) 


4.32(23) 


step width 


39.253(45) 


47.781(73) 


49.881(96) 


44.01(18) 


4.31(24) 


running width 


39.251(49) 


47.781(74) 


49.898(95) 


44.48(22) 


10.83(28) 


complex mass 


39.221(45) 


47.730(73) 


49.770(91) 


43.97(18) 


4.31(23) 



Table 4: Cross sections for the processes 77 — > e z/ e z/ M /i + and 77 — > e z/ c z/ M/ u + 7 for various 
CM energies and various width schemes without convolution over the photon spectrum. 



show that a naive signal definition is a bad concept for 77 — > WW — > 4/, since deviations 
from the full process 77 — > 4/ even reach 5-10% in the TeV range. This is in contrast to 
the situation at e + e~ colliders where the naive W-pair signal (defined in 't Hooft-Feynman 
gauge) was a reasonable approximation (see, e.g., Ref. [35]). 

The failure of the naive W-pair signal definition for 77 collisions was also pointed out 
in Refs. [6,7] before. In Ref. [7] an "improved narrow-width approximation" was presented 
which provides another variant for a gauge-invariant W-pair signal definition. It is based 
on the factorization of production and decay matrix elements, while retaining W-spin 
correlations. 

6.2.5 Comparison of schemes for introducing finite gauge-boson widths 

In this section we compare the different implementations of gauge-boson widths de- 
scribed in Section 2.4 numerically. As explained in Section 2.4, the complex-mass scheme 
is the only scheme that yields gauge-invariant results in general, but for the process classes 
77 — > 4/(7) the fixed- width approach (in the non-linear gauge) also yields amplitudes 
that respect Ward identities and gauge cancellations. Table 4 lists the cross sections for 
the processes 77 — > e~z/ c z/ M/ u + and 77 — > e~z/ e z/ At /i + 7 obtained with the fixed W width, 
the step-width, the running-width, and with the complex-mass scheme. The results of 
all four schemes for the process 77 — > e~u c u^ii + agree within the expected accuracy of 
0(T-w/M-w) up to energies in the TeV range. However, for 77 — > e~z/ e ^Ai + 7 the running- 
width scheme yields totally wrong results for several TeV, while the other schemes are still 
in good agreement. Although the gauge-invariance-breaking effects in the running-width 



27 




Figure 8: Invariant-mass distribution of the W + boson reconstructed from the ud quark 
pair (l.h.s.) as well as its pro duct ion- angle distribution (r.h.s.) in the reaction 77 — > 
e~u e nd at = 500 GeV with and without convolution over the photon spectrum. 



scheme are formally of 0(Tw/My/), they are enhanced by spoiling gauge cancellations, 
thereby ruining the reliability of the prediction completely 

For the semileptonic 77 — > 4/ process it was already observed in Ref. [7] that the 
cross section does not vary significantly if the fixed-width, the running-width, or a so- 
called "fudge-factor" scheme is used for introducing finite widths. 

6.3 Distributions 

6. 3. 1 Distributions for reconstructed W bosons 

In Figure 8 we show the invariant-mass and angular distributions of the intermediate 
W + boson for the process 77 — > e~z/ e ud at y/s = 500 GeV. The momentum of the W + 
boson is reconstructed from the outgoing quark pair in the decay W + — > ud. Figure 8 
also illustrates the effect of the convolution over the photon spectrum. 

The resonance in the invariant-mass distribution (l.h.s. of Figure 8) has the typical 
Breit-Wigner shape and can be used to determine the W-boson mass and width at a 77 
collider. Moreover, owing to its large cross section, the W reconstruction in this reaction 
seems to be a promising possibility for detector calibration at a 77 collider. Similarly to 
the integrated cross sections discussed in the previous sections, the convolution qualita- 
tively rescales the distribution by roughly a factor 4. 

The r.h.s. of Figure 8 shows the distribution in the angle # u g between the W + boson 
and the beam axis. Since the incoming 77 state is symmetric w.r.t. interchange of the two 
photons, the angular distribution is symmetric in the production angle 8 U ^. W bosons are 
predominantly produced in forward or backward direction owing to diagrams with t- and 
li-channel exchange of W bosons. For the process 77 — > WW with on-shell W bosons, 
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the forward and backward peaks are integrable and lead to a constant cross section in 
the high-energy limit. As already pointed out in Section 6.2.2, the angular cuts (6.2) 
restrict the available phase space of the intermediate W bosons and lead to a reduction 
of the forward and backward peaks for sufficiently high energies. Note that the reduction 
induced by the convolution over the photon spectrum is not uniform, but tends to flatten 
the shape of the angular distribution slightly. This is mainly due to the reduced CM 
energy in the photon spectrum, leading to a less pronounced peaking behaviour in the 
forward and backward directions. 

6.3.2 Energy and production- angle distributions of fermions 

In Figure 9 we show the energy and angular distributions of the outgoing fermions e~, 
u, and d in the reaction 77 — > e~z/ c ud at y/s = 500 GeV with and without convolution 
over the photon spectrum. 

For monochromatic, unpolarized incoming 7 beams (i.e. without convolution over the 
photon spectrum), the energy distributions (l.h.s. of Figure 9) of the fermions e~, u, and 
d almost coincide and are maximal at their largest and smallest kinematical limits. These 
regions are dominated by the situations where the respective W boson emits the considered 
fermion parallel or anti-parallel to its direction of flight. The convolution over the photon 
spectrum changes the shapes of the energy distributions considerably. Since the photon 
spectrum falls off rapidly for energies above 80% of the incoming electron energy, energies 
of the final-state fermions larger than 200 GeV become practically impossible. For fermion 
energies below 200 GeV the shapes of the distributions of the outgoing fermions e~ and 
u look rather different from the one for the anti-fermion d. This effect is due to the 
effective 7 beam polarization in the photon spectrum; for unpolarized 7 beams the energy 
distributions would look almost identical. In detail, the effective polarization of the 77 
system is mainly (A1A2) = (++), leading predominantly to W + W~ production with 
effective helicities (++) (see, e.g., Refs. [7,21]). Owing to helicity conservation, however, 
W bosons with helicity +1 cannot decay into fermion-anti-fermion pairs with a fermion 
(which must have helicity — \) parallel to the flight direction of the W boson. Thus, 
much more anti-fermions (which have helicity +|) than fermions follow the directions of 
the decaying W bosons, which qualitatively explains the reduction (enhancement) of the 
fermion (anti-fermion) energy distributions at the upper kinematical energy limit. The 
above arguments are nicely illustrated in Ref. [7], where the fermion energy distributions 
are shown for fully polarized, monochromatic photon beams. 

The r.h.s. of Figure 9 shows the distributions in the angles Of of the ( ant i-) fermions 
/ = e~,u,d to the beam axis. Because of the symmetry of the incoming 77 state w.r.t. 
interchange of the two photons, the angular distribution is symmetric in Of. The forward 
and backward peaks originate from two sources. The by far largest contribution to the 
differential cross section comes from signal diagrams and thus from configurations where 
the W bosons as well as the decay fermions are nearly parallel to the beam. The second 
source, which is widely suppressed by the applied cuts, is related to collinear singularities 
of background diagrams where an incoming photon splits into an fermion-anti-fermion 
pair //, with the fermion or anti-fermion directly going into the final state. If the phase 
space of the outgoing ( ant i-) fermion is not restricted by cuts, such collinear or mass 
singularities lead to logarithms of the form \n(s/m 2 f), where m/ is the fermion mass. Since 
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Figure 9: Energy (l.h.s.) and production-angle (r.h.s.) distributions of the outgoing 
fermions e~, u, and d in the process 77 — ► e~z/ e ud at y/s = 500 GeV with and without 
convolution over the photon spectrum. 
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Figure 10: Invariant- mass distribution of the four-quark final state for the process 77 
udsc at y/s^ = 260 GeV including convolution over the photon spectrum. 



our calculation is done for massless fermions, the collinear singularities must be excluded 
by phase-space cuts and the fermion mass in In^s/mj) is replaced by the corresponding 
cut parameter. 

The photon spectrum reduces the differential cross section over the whole range and 
again flattens the angular distributions, especially in the cases of outgoing fermions. The 
significant difference between the outgoing fermions and anti-fermions is again due to the 
effective 7 polarization in the photon spectrum. As explained above, more anti-fermions 
than fermions follow the flight directions of the W bosons, which are mainly produced in 
the forward and backward directions. This is the reason why the 8 e - and 8 U distributions 
are flattened, while the peaking behaviour in the 8^ distribution is more pronounced after 
the convolution over the photon spectrum. 

6.3.3 Higgs-boson resonance 

In Figure 10 we show the invariant-mass distribution of the Higgs boson for the process 
77 — > H — > WW — > udsc for a Higgs mass of Mh = 170 GeV. The CM energy of the 
electron beams is chosen to be y/s^ = 260 GeV which maximizes the 77 luminosity in the 
region y/s^ ~ Mr. The invariant mass M u ^ s5 of the Higgs boson is reconstructed from 
its decay products which are the four outgoing quarks. This means that M ud j sE is equal to 
the photonic CM energy, M U( j sE = ^/s 77 . Thus, the shape of the distribution depends on 
the form of the photon spectrum very strongly. The effective 77H coupling is set to the 
SM value (4.3). For comparison the situation without Higgs resonance is also included in 
Figure 10, illustrating the significance of the Higgs signal. The different peak heights in 
the two plots simply result from different bin sizes. 
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Figure 11: Combined cross section for semi-leptonic final states as a function of the ATGC 
Ak 7 and A 7 (l.h.s.) and lcr contours (r.h.s.) in the (Ak 7 , A 7 ) plane at ^/s^ = 500 GeV 
including the convolution over the photon spectrum. 



6.4 Anomalous couplings 

In this section we study the impact of possible anomalous gauge-boson couplings on 
CC cross sections of the process class 77 — > 4f. In order to estimate the full sensitivity of 
a future 77 collider, such as the 77 option at TESLA, on anomalous couplings, in addition 
differential distributions and realistic event selections should be taken into account. Such 
a study goes beyond the scope of this paper, but our Monte Carlo generator can serve as 
a tool in this task. 

We consider only semi-leptonic final states, since these have the cleanest experimental 
signal. The cross section for semi-leptonic final states is obtained from the sum over all 
reactions 77 — > l~viqq', with q = u, c and I = e, /1, r, and their corresponding charge- 
conjugated processes 77 — > vil + q'q. The results are shown in Figure 11 for ATGC and in 
Figure 12 for AQGC. In the left plot of Figure 11 and the upper plot of Figure 12 we show 
the cross section as a function of the anomalous coupling constant normalized to the SM 
cross section. As can be seen in the insert of Figure 11, the minimum in the Ak 7 curve is 
shifted to negative values which is caused by contributions to the cross section that are 
linear in Ak 7 . These contributions result from the interference between matrix elements 
linear in the ATGC Ak 7 with the SM amplitude. On the other hand, the interferences 
for the ATGC A 7 are small. In the case of AQGC, such interferences are relatively large 
for a c . 

In order to examine the sensitivity of a linear collider to anomalous couplings, we 
consider a 77 collider with an integrated luminosity of C = 100 fb~ and a CM energy of 
^ = 500 GeV [1]. We define 

X 2 = {N{at) ~ N)2 with N = a SM C, N(at) = *(<*)£, (6.4) 
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Figure 12: Combined cross section for semi-leptonic final states as a function of the 
AQGC ao, a c , and (upper plot) and la contours (l.h.s. projection, r.h.s. section) in the 
coordinate planes at y/s^ = 500 GeV including the convolution over the photon spectrum. 
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where N is the expected number of events in the SM and iV(aj) the number of events in 
the SM extended by the non-standard couplings. In Figures 11 and 12 the la contours 
corresponding to x 2 = 1 are shown. Note that the la contour can result from iV(aj) > iV 
and N(di) < N. While x 2 = 1 with iV(aj) > iV is always possible for sufficiently large 
anomalous couplings, = 1 with N{ai) < N requires large interference effects of matrix 
elements with anomalous couplings. In our case, both branches of the la contours are 
realized. In Figure 11 the plot on the r.h.s. shows the la contours in the (A« 7 , A 7 ) plane. 
Since the cross section is a polynomial up to fourth power in the ATGC, the contours 
are not of elliptic form. The allowed region lies between the two contours that are rather 
close to each other so that they cannot be distinguished in the insert which shows the 
contours on a larger scale. Note that in the limit of large luminosity the contour in the 
insert of the r.h.s. of Figure 11 does not shrink to a point, but reduces to a line in the 
(A/c 7 , A 7 ) plane on which <x anom — °sm- I n order to resolve this correlation between Ak 7 
and A 7 , anomalous effects on distributions should be considered, or other constraints from 
the e + e~ or e~7 modes should be included. 

In case of AQGC the cross section is at most quadratic in the AQGC, and the x 2 — 1 
surface consists of two ellipsoids in the (a ,a c ,a ) space. The existence of two branches 
is again due to large interferences of anomalous contributions. In the lower left plot of 
Figure 12 we show the projections of the outer ellipsoid into the coordinate planes of two 
AQGC (where the third AQGC is zero). In the lower right plot the sections of both 
ellipsoids with these planes are given. Since the centre of the ellipsoids is shifted in the a c 
and ao directions, the terms in a(cii) linear in these couplings are significant; they result 
from interferences of the diagram with the AQGC with the SM amplitude. Interferences 
that are proportional to a turn out to be small. From Eq. (3.14) it is obvious that there 
are no aoao and a c ao terms in a(a,i). Consequently, the projection into and the section 
with the (ao, a c ) plane coincide. On the other hand, the two other projections and sections 
differ, signalling that the a c ao term in cr(aj) is significant. 

The allowed la region (x 2 < 1) in the (ao, a c , ao) space is the shell at the boundary of 
the shown ellipsoid. Similar to the observation made above for the ATGC, the size of the 
ellipsoid does not shrink for larger luminosity, only the thickness of the shell will decrease. 
This means that the size of the projections shown in the lower left plot of Figure 12 will not 
reduce for larger luminosity. Thus, using only information on an integrated cross section 
(for a fixed energy) could not improve the bounds on AQGC w.r.t. the ones resulting 
from e + e~ — > WW7 — > 4/7 [32]. However, the thinness of the shell of the ellipsoid, as 
illustrated in the lower right plot of Figure 12, shows that the bounds can be drastically 
tightened if the correlation between the three AQGC is resolved. Differential distributions 
will certainly provide this information, so that a 77 collider should be able to constrain 
AQGC by an order of magnitude better than an e + e" collider operating at comparable 
energy. 

7 Summary 

A Monte Carlo generator, which is based on multi-channel techniques with adaptive 
self-optimization, has been constructed for lowest-order predictions for the processes 77 — > 
4/ and 77 — > 4/7. Fermions are treated in the massless approximation. Anomalous triple 
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and quartic gauge-boson couplings, as well as an effective 77H coupling, are included in 
the transition matrix elements for 77 — > 4f. For 77 — > 4f compact results for the helicity 
amplitudes are presented in terms of Weyl-van-der-Waerden spinor products. 
Using this generator, we have presented a variety of numerical results: 

• A representative list of integrated cross sections for the processes 77 — > 4/ and 
77 — > 4/7 is compared to results obtained with the Whizard/Madgraph package. 
We find agreement between both Monte Carlo programs. 

• The dependence of some 77 — > 4/ and 77 — > 4/7 cross sections on the CM energy 
is shown. In this context, the influence of a realistic photon beam spectrum and the 
size of subcontributions originating from charged-current, neutral-current, or strong 
interactions are investigated. 

• The complete lowest-order cross section for 77 — > WW — > 4/ is compared with 
the corresponding double-pole approximation, revealing an accuracy of the latter of 
1—3%. While this deviation is of the naively expected order of r w /M w , the gauge- 
variant subset of diagrams with two resonant W bosons, the so-called "CC03" signal 
diagrams, yields differences of 5—10% to the full calculation. The deviations of the 
naive signal grow with increasing energy, whereas the quality of the double-pole 
approximation remains stable. 

• Different schemes for introducing a finite width for gauge bosons are analyzed, 
and good agreement between the gauge-invariant complex-mass scheme and the 
fixed-width scheme is found. However, the results also reveal problems (at least in 
77 — > 4/7) when using the running width. 

• A few differential distributions are discussed for 77 — > 4/ processes that proceed via 
charged-current interactions, in particular comprising distributions in the invariant 
mass and in the production angle of reconstructed W bosons and in the invariant 
mass of a resonant Higgs boson. Moreover, it is shown that the convolution over 
the photon spectrum significantly distorts energy and angular distributions of the 
produced fermions due to an effective photon polarization. 

• Finally, we examine the effects of anomalous triple and quartic gauge-boson cou- 
plings on 77 — > 4/ cross sections. Since contributions of anomalous couplings to 
cross sections can cancel in specific configurations, it is necessary to take into ac- 
count results from other observables (such as differential distributions) or from other 
experiments (such as e + e~ or e~7 collisions) in order to constrain individual anoma- 
lous couplings. Our results suggest that an analysis of the processes 77 — > 4/ can 
constrain anomalous 77WW couplings about an order of magnitude better than 
studying e + e~ — > 4/7. 

Apart from studying the processes 77 — > 4/ and 77 — > 4/7 in lowest order the 
described event generator is a first step towards a precision calculation of the process 
77 — > WW — > 4/ which has to include radiative corrections to the dominating W-pair 
production process. In particular, the results for the radiative processes 77 — > 4/7 can 
serve as a basic building block in this task. In a next step, we will complete the calculation 
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of electroweak corrections in the double-pole approximation similar to the approach of 
RacoonWW for e + e~ annihilation. 
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